Quantum Statistics of Interacting Dimer Spin Systems 
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The compound TICUCI3 represents a model system of dimerized quantum spins with strong in- 
terdimer interactions. We investigate the triplet dispersion as a function of temperature by inelastic 
neutron scattering experiments on single crystals. By comparison with a number of theoretical ap- 
proaches we demonstrate that the description of Troyer, Tsunetsugu, and Wiirtz [Phys. Rev. B 50, 
13515 (1994)] provides an appropriate quantum statistical model for dimer spin systems at finite 
temperatures, where many-body correlations become particularly important. 
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The temperature is without doubt the key control pa- 
rameter in solid-state physics, for both historical and 
technical reasons. Historically, the study of thermal 
fluctuations created the discipline of statistical mechan- 
ics. As one technical example, the liquefaction of he- 
lium in 1908 enabled the discovery of superconductivity 
in 1911 [l|. While control parameters such as the mag- 
netic field or applied pressure may seem rather abstract 
by comparison, in certain classes of system they have 
a profound effect on the quantum mechanical fluctua- 
tions of the dominant degrees of freedom. In this respect, 
dimer spin systems such as BaCuSi2 0g 000, TlCuCl 3 
00 011 El and Cs3Cr 2 Brg 01 have recently attracted 
considerable attention as they exhibit both field- and 
(in TICUCI3) pressure-induced magnetic quantum phase 
transitions. The ordered phases have also been described 
as a novel state resulting from Bose-Einstein condensa- 
tion of magnons 0, . 

To date the primary experimental and theoretical 
studies of these systems have focused on elucidating 
the parameter-dependence of the ground and excited 
states at "T = K" 03. However, 
the temperature-dependence of the spin-spin correlation 
function S(Q,u>) in such dimer-based compounds is also 
of fundamental interest. At finite temperatures, thermal 
fluctuations populating the excited triplet states above 
the energy gap, A w 7.5 K, in the spin excitation spec- 
trum of TICUCI3 become significant and compete with 
the intrinsic quantum fluctuations. The absence of a 
classical ordering transition in zero field and at ambient 
pressure, and the strong interdimer interactions, make 
TICUCI3 an ideal candidate for an experimental investi- 
gation of the quantum statistical description for such a 
spin system. 



In TICuCFj, the antiferromagnetic (AF) dimer unit is 
formed by the two 5 = 1/2 magnetic moments in a pair 
of Cu 2+ ions, and has a singlet ground state (total spin 
5 = 0) with triplet (5 = 1) excited states which may be 
considered as bosonic quasiparticles. Only in the high- 
temperature limit, relative to the intradimer exchange 
constant J rj 56 K ^J, is the system better de- 
scribed in terms of 5 = 1/2 fermions whose interactions 
are weak compared to the thermal fluctuations, and AF 
dimer correlations are strongly suppressed. In the inter- 
mediate regime, A < T < J, a singlet-triplet descrip- 
tion remains valid but the number of excited triplets 
becomes significant. Here the effect of interparticle in- 
teractions, particularly the intrinsic hard-core repulsion 
on each dimer bond, defines a challenging experimental 
and theoretical scenario. Temperature-dependent prop- 
erties of dimer spin systems have to date been measured 
only in compounds of weakly interacting dimers with a 
spin gap large compared to the triplet bandwidth, includ- 
ing BaCuSi 2 6 0, KCuCl 3 H3, and Cu(N0 3 )2-2.5D 2 
|l4j . where many-body correlation effects are less pro- 
nounced. The description of these data was restricted to 
RPA theory or phenomenological modeling. 

In this Letter we present the results of a recent inelas- 
tic neutron scattering (INS) study of the spin dynam- 
ics in TICUCI3. The primary results include a compre- 
hensive experimental survey of the triplet excitations up 
to temperatures approaching the intradimer interaction 
strength, A <C T max ~ 0.7 J, and their theoretical model- 
ing by quantum statistical considerations applied to the 
known zero-temperature states (i.e. with no additional 
fitting parameters). We characterize the evolution of the 
spin gap, A(T), which is of particular relevance to the 
quantum critical behavior of the compound, and shows a 
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pronounced increase with temperature. We demonstrate 
that the data support strongly the Troyer-Tsunetsugu- 
Wiirtz (TTW) Ansatz for the quantum statistics of gen- 
eral dimer spin systems (lij . Our results may further 
be compared with those obtained for weakly interacting 
dimer compounds (above) and are closely related to ex- 
tensive investigations of another class of gapped quantum 
magnets, S = 1 Haldane chains, which have been studied 
in model compounds such as CSMCI3 

Thermal renormalization of the dispersion, spectral 
weight and damping of the triplet modes in KCUCI3, 
which shares the crystal structure and exchange coupling 
scheme of TICUCI3, are adequately described within a 
self-consistent RPA theory including higher-order cor- 
rections ^jj ■ The considerably stronger interdimer inter- 
actions make this approach less successful for TICUCI3. 
We have therefore applied an extension of the bond- 
operator description of the compound 0] to finite tem- 
peratures, a procedure requiring [l8j | the inclusion of 
thermal population factors n(E a (q), [3) in the mean- 
field (MF) equations defining the singlet density s and 
the triplet chemical potential fj,; here E denotes the 
magnon energy, q its wave vector, (3 = {k^T)~ l , and 
a = — 1,0, +1 represents the S z — — 1,0, +1 triplet 
modes. We assume explicitly that there are no addi- 
tional factors which may change with temperature, such 
as renormalization of the intra- and interdimer interac- 
tion parameters due to phonons 

The singlet state and each of the triplets on a single 
dimer obey an exclusion constraint and are described ap- 
propriately as hard-core bosons. Their statistical proper- 
ties have no straightforward extension to reciprocal space 
and no exact expression exists for a triplet thermal oc- 
cupation function. However, the effects of the constraint 
were introduced by Troyer and coworkers [l5| in the form 
of a statistical reweighting of the number of free-boson 
states of the iV-dimer system to account for their com- 
plete local exclusion. The reweighting constitutes a sup- 
pression of state availability in all magnon occupation 
sectors M, being most significant at large M, and re- 
turns the effective single-dimer free energy 
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with z a (0) — jfJ2q e ~ >3Ea ^ the partition function for 
triplet mode a. The effective occupation function for this 
statistical Ansatz (denoted TTW-MF) is then 



n(E a (q),(3) 



e -/3E a (q) 

1 + Ea *«(/?)' 



(2) 



Equation (2) reproduces both the low- and high- 
temperature limits expected for the constrained boson 
system ^ij, namely n ~ e~^ Ea ^ at low T (bosonic, M 
small) and n ~ 1/4 at high T (equipartition, M large), 
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FIG. 1: (a) Finite-T triplet dispersion for TIC11CI3 from sta- 
tistical models described in the text. A denotes the mini- 
mum singlet-triplet gap, Q~q+T the momentum, where r is 
a reciprocal-lattice vector, and open circles the Q values pre- 
sented in Fig. |21 (b) Temperature-dependence of A(T) from 
TTW-MF, Bose-MF, and RPA approaches. 



and we will test below the accuracy of exclusion imple- 
mentation by statistical reweighting for the regime be- 
tween these limits. For comparison we will show also 
the results (denoted Bose-MF) obtained by enforcing the 
constraint at the Holstcin-Primakoff level and con- 
sidering a conventional bosonic magnon occupation fac- 
tor. Fig. Ufa) shows theoretical predictions for the ther- 
mal renormalization of the triplet dispersion: both ap- 
proaches display a narrowing of the magnon band which 
arises because interdimer triplet hopping is obstructed by 
the exclusion condition when the concentration (M/N) 
of thermally excited triplets becomes significant. 

The gap A at the dispersion minimum, which becomes 
the AF zone center in the magnetically ordered high- 
field and -pressure phases, deserves special attention. 
Fig. Gib) shows its thermal evolution as obtained from 
the above models, and from the RPA treatment applied 
to KCUCI3 in Ref. . The Bose-MF approach shows a 
very strong renormalization at higher temperatures due 
to the possibility of unlimited boson occupation. This 
effect is considerably weaker in the TTW-MF calcula- 
tions, where it is limited by the exclusion encoded in the 
statistics. These approaches coincide at low tempera- 
tures, where the triplet occupation is too dilute for their 
hard-core nature to be relevant. By contrast, in the RPA 
treatment the onset of a perceptible bandwidth narrow- 
ing requires much higher temperatures. 

These theoretical results clearly motivate INS mea- 
surements of the triplet dispersion relation as a function 
of temperature at zero and finite field. We have con- 
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FIG. 2: (a) INS spectra for TlCuCla measured on IN14 at 
Q=(0 4 0) r.l.u. Solid lines are Lorentzian-based least- 
squares fits including full 4D-resolution convolution (peak 
asymmetry) and triangles indicate the physical excitation en- 
ergy, (b) and (c) INS intensity contour plots measured at 
1.6 K and 39.8 K, compared with theoretical curves by the 
Bose-MF (dashed) and TTW-MF (solid lines) descriptions. 



ducted INS experiments on single crystals of TICUCI3 at 
H = 0. These were performed on the triple-axis spec- 
trometers TASP (SINQ), IN14 (ILL), and 2T (LLB), 
working in constant-fc/ mode with a focusing PG ana- 
lyzer/monochromator. While 2T is operated with ther- 
mal neutrons of kf = 2.662 A -1 (14.7 meV) and 1.970 
A" 1 (8.0 meV), TASP and IN14 use moderated (cold) 
neutrons with a characteristic spectrum shifted to lower 
energies, kf = 1.506 A" 1 (4.7 meV). A cooled beryl- 
lium or PG filter is positioned between the analyzer and 
the sample, contained in a standard cryostat providing 
T > 1.5 K. The (horizontal) collimation adopted on the 
cold instruments was 60'-open-open and 40'-open-open, 
yielding an energy resolution of 0.2 meV (full width at 
half maximum height) in both cases. Corresponding val- 
ues for the thermal instrument without additional colli- 
mation were 1.2 meV and 0.5 meV. 

Typical INS spectra for the triplet modes in TICUCI3 
measured for a number of temperatures at the band min- 
imum (Q= (Q 4 0) r.l.u.) are presented in Fig.^a). Up- 
ward renormalization of the excitation energy, the onset 
of damping, and reduction in the integrated intensity are 
observed. Figs. EJb) and (c) contrast the INS intensi- 
ties measured at 1.6 and 39.8 K along (Qh 0). The 
complete experimental survey, including representative 
Q values spanning the entire bandwidth, is presented in 
Fig. Eland compared with the model calculations. For in- 
strumental reasons the spin gap A must be measured in 
several experimental configurations to cover the full tem- 
perature range. The point Q = (1.35 0) r.l.u. is of par- 




TTW-MF 
Bose-MF 
2T (therm.) 
TASP (cold) 
IN14 (cold) 



10 15 20 25 30 35 40 
Temperature [K] 

FIG. 3: Temperature-dependence of the excitation energies 
at (I) Q=(0 1) (closed black circles and closed red circles) 
and Q=(0 4 0) (open black circles), (II) Q=(0.25 -1.5), (III) 
Q=(0.5 -2), (IV) Q=(1.35 0), and (V) Q=(1.25 0) (all 
in r.l.u.). Data from different instruments as indicated. Solid 
lines are obtained from TTW-MF and Bose-MF calculations. 



ticular interest: here the excitation energy corresponds 
approximately to J, and indeed it is nearly temperature- 
independent. The downward renormalization of the dis- 
persion maximum is evident. The continuity between 
data sets from different instruments and Q values under- 
lines the accuracy of the experimental procedure, includ- 
ing the data analysis (below). 

With regard to the calculated fits, the TTW-MF 
model is clearly favored by the experimental observa- 
tions in the range 20 K < T < 40 K (Figs. Etc) and EJ 
where the thermal renormalization changes less rapidly 
due to the constrained magnon population. The Bose- 
MF approach fails markedly in this respect. For clarity 
the RPA predictions are not shown in Fig. El and we 
state only that these fail in particular to reproduce the 
low-temperature behavior [cf. Fig.^b)]. The agreement 
between the TTW-MF result and the experimental data 
is excellent over the full temperature range of the investi- 
gation. We stress that no free parameters are introduced 
with respect to the low-temperature dispersion |5ll7|. 

We turn to a discussion of the excitation line width. 
For modes of infinite lifetime, S(Q,uj) is a rj-function in 
energy (E = hcu), which within instrumental resolution 
is the case for the triplet modes measured in TICUCI3 
at T = 1.6 K. However, at elevated temperatures a fi- 
nite intrinsic line width (damping) is observed, and data 
analysis requires knowledge of the spectral function. Dif- 
ferent 0(n) rotor models give a Lorentzian line shape in 
one dimension (ID) and at T < A [i^, but neither limit 
is appropriate here. By contrast, the damped harmonic 
oscillator (DHO) line shape used frequently for phonons 
[2l| is identical to a double Lorentzian if the real ex- 
citation energy is defined as -Edho('Z) 2 = E(q) 2 + T 2 , 
and thus depends on the damping T (half Lorentzian 
width at half maximum height). The DHO definition 
has been used to fit the finite-T triplet excitations in 
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FIG. 4: Temperature-dependence of line width IT for selected 
reciprocal-lattice points. Solid lines denote fits to the data 
as shown in the legend (A=3(l), C=2.7(5), and m=2.04(7)). 
The low-T offset of 0.05 meV is attributed to uncertainty in 
the resolution deconvolution. 



Haldane chains 0], and also for spin waves |2^ . 

In the Bose-MF and TTW-MF approaches the energy 
is not renormalized by line shape effects, consistent with 
the Lorentzian assumption. While a line width T may be 
computed from the phase space available for magnon de- 
cay and scattering, these considerations provide no proof 
of either a Lorentzian or a DHO form for triplet exci- 
tations in 3D, and nor is one available in the literature. 
We proceed on the basis of experiment: the INS mea- 
surements show that a Lorentzian yields an accurate de- 
scription of the data. The spectra [Fig. |2Ia)] were fitted 
with this line shape and the real excitation energy de- 
fined as E(q). We suggest that this definition provides 
the correct description of the spectral function at finite 
temperatures in dimer spin systems. This conclusion is 
reinforced by the temperature-dependence of the disper- 
sion maxima, which show a downward renormalization 
with increasing T. This behavior is consistent with the 
Lorentzian model, whereas -Edho(<z) describes the max- 
ima rising due to the strong increase of T with tempera- 
ture. Qualitatively, this result may be traced once again 
to the hard-core nature of the magnon excitations. 

Fig. 01 shows the damping of the triplet modes, which 
is clearly activated. This behavior is expected from any 
model based on enhanced triplet-triplet scattering as the 
origin of a finite quasiparticlc lifetime. We are unable to 
distinguish between fits using a variable activation energy 
or a variable power-law prefactor to an exponential with 
a constant activation energy of J |3 (Fig- LH ■ However, 
the data is definitely incompatible with an activation en- 
ergy of 2A(T), as would be expected from magnon decay 
terms present at third order in the triplet hopping Hamil- 
tonian, and these processes may be excluded. 

With the present investigation we would hope to ini- 
tiate further studies. Detailed theoretical predictions for 
the finite-temperature spectral function of gapped exci- 
tations in 3D dimer spin systems are required. Experi- 



mental studies at magnetic fields around the critical field 
H c and at finite temperatures, i.e. in the quantum critical 
regime, are further expected to reveal new physics. 

In summary, we have investigated the thermal renor- 
malization of triplet excitations in the strongly interact- 
ing dimer spin system TICUCI3 by means of INS. We find 
excellent agreement between the TTW-MF Ansatz and 
the complete data set, which was taken up to T « 0.7 J. 
The Bose-MF approach is in comparable agreement for 
T < 20 K and may therefore be considered as a suitable 
low-T approximation. However, at higher temperatures, 
where many dimer triplets are excited, their hard-core 
nature becomes manifest and an accurate quantum sta- 
tistical description requires the TTW Ansatz. 
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